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ABSTRACT 



According to our present understanding, long gamma-ray bursts (GRBs) originate from the collapse of massive stars, while short 
bursts are caused by to the coalescence of compact stellar objects. Because the afterglow evolution is determined by the circumburst 
density profile, n(r), traversed by the fireball, it can be used to distinguish between a constant density medium, n(r) = const., and 
a free stellar wind, n(r) cc r" 1 . Our goal is to derive the most probable circumburst density profile for a large number of Swift- 
detected bursts using well-sampled afterglow light curves in the optical and X-ray bands. We combined all publicly available optical 
and SwiftfK-my afterglow data from June 2005 to September 2009 to find the best-sampled late-time afterglow light curves. After 
applying several selection criteria, our final sample consists of 27 bursts, including one short burst. The afterglow evolution was then 
studied within the framework of the fireball model. We find that the majority (18) of the 27 afterglow light curves are compatible with 
a constant density medium (ISM case). Only 6 of the 27 afterglows show evidence of a wind profile at late times. In particular, we 
set upper limits on the wind termination-shock radius, R T , for GRB fireballs that are propagating into an ISM profile and lower limits 
on R T for those that were found to propagate through a wind medium. Observational evidence for ISM profiles dominates in GRB 
afterglow studies, implying that most GRB progenitors might have relatively small wind termination-shock radii. A smaller group of 
progenitors, however, seems to be characterised by significantly more extended wind regions. 

Key words, gamma-ray burst: general, ISM: structure, Radiation mechanisms: non-thermal 



1. Introduction 

Starting with the discovery of the gamma-ray b urst- supernova 
(GRB-SN) association GRB 980425/SN 1998bw fcalama et alj 
1998; llwamotoetai]ll998t ISollerman et all l20od) . there is bv 
now convincing evidence that the progenitors of long GRBs 
are massive stars exp l oding as type Ic SNe (for a review 
see IWooslev & Blo om 2006). Within this picture, the opti- 
cal light observed after a long GRB is the superposition of 
the afterglow light, a supernova component, and light from 
the underlying host galaxy (plus potential additional radiation 
components at very early times, which we will not consider 
here). Phenomenologically, this immediately unveils two ob- 
serving strategies to reveal a massive-star origin of a GRB: 
(z) via the detection of a late -time SN bump (method i) 
in the optical light curve (e.g., Reichartl Il999t iGalama et alj 
2000; iDado et alll200l IZeh et al] 12001 and (ii) by the spec- 
troscopic confirmation (method ii) of associate d supernova 
light (the best case so far being GRB 03 0329: Hiort h et al 



the G RB progenitor (e.g.. iMirabal et alj l2003t ISchaefer et al.1 
20031; iKlose et alJl2004l IStarhng et alj|2005h iBerger et alj|2006; 
Fox et al.l 12008: Cas tro-Tirado et a02010l) . Some authors, how- 



2003 



Kawa bata etaDl2003t iMatheson etai]l2003l: [Stanek et al 



2003). In addition to both observing strategies, there are two 
further methods by which a massive-star origin can be re- 
vealed. Some afterglow spectra showed blue-shifted absorption 
line systems (method Hi), which can be understood as signa- 
tures from the expanding pre-explosion wind escaping from 



ever, notice that several properties of the putative blue-shifted 
absorption line systems disagree with the expectations from 
Wolf-Rayet ( WR) winds, e.g., line widths, ionization levels and 
metallicities dChen et alj|2007t iProchaska et al.ll2007t : iFox et al] 
2008). Finally (method iv), the circumburst medium determines 
the spectral and temporal evolution of the afterglow, allowing us 
to discern between a constant- density medium , n(r) _ const., 
and a wind medium, n(r) oc r~ 2 (ISari et alj [T998; Chevalier & Li 
1999. l2000h. This meth od was successfull y applied in, e.g., 
Istarling et all (120091) and lCurran et al.l (120101) . 

Naturally, these various approaches have their observa- 
tional advantages and disadvantages. While method (ii) can 
provide the strongest observational evidence for a massive-star 
origin of the GRB under consideration, it can only be ap- 
plied to the nearest and hence brightest events up to a red- 
shift of about 0. 5. Even 13 years after the first di s cover y 
of an afterglow (ICosta etalJ 1 1 9971: Ivan Paradiis et all Il997l) . 
i.e., after more than 500 GRBs with detected (X-ray, opti- 
cal, radio) afterglow lighfl secure evidence for a spectroscop- 
ically associated SN was only reported f or rough ly 1% of all 
event s (GRBs 980425:lGalama et al.lll99l 03 0329: lHiorfh et al. " 
120031: iKawabata et al.ll2003b IMatheson et al.ll2003t IStanek et al 
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2003 1 : 031203: iMalesani et alJ l200fl 060218: iFerrero et al] 
2006; Mi rabal et alJ 120061: iModiaz et alJ 120061: iPian et all 
2006; ISollerman et alJ 120061; 081007: iDella Valle et al.l 12008; 
1 003 1 6D: IChornock et alJl2010t IStarling et al.ll2010h . Contrary 
to this approach, method (/) can basically reveal a SN compo- 
nent up to a redshift of about 1 CZeh et ai]|2004l) . assuming a 
non-extinguished SN 1998bw as a template (the m ost distant SN 
bump was found for GRB 000911 at z = 1.06: iMasetti et all 
12005b - At notably higher redshifts a GRB-SN becomes too faint 
to be discovered even with 8m-class op tical telescopes because 
of line blanketing (e.g., Filippenko 1993). 

In principle, methods (z'z'z) and (z'v) do not have redshift con- 
straints, because both rely on the observation of the afterglow 
and not on the (expected) SN component. In addition, method 
(z'v) splits into different approaches and basically works for the 
optical and X-ray band in the same way. Moreover, it can even 
be applied to the most distant GRBs, which are already affected 
by Lyman dropout in the optical bands. Its main disadvantage is 
that it usually requires substantial observational efforts. In par- 
ticular, it relies on the measurement of the light-curve evolution, 
including the determination of the spectral energy distribution 
(SED). 

In this paper, applying method (iv), we use late-time data of 
afterglows detected by Swift to tackle the question of the pre- 
ferred density profiles in a statistical sense. Our goal is to com- 
bine all publicly available optical data with Swift/XRT data in or- 
der to determine the corresponding circumburst density profile, 
n(r) oc r~ k . Qualitatively, this splits into either a constant-density 
medium (hereafter referred to as the interstellar medium or ISM) 
(k = 0) or a free wind profile (k = 2). While the complex wind 
history of an evolved massive star might p roduce density pro- 
files different from the ideal case (k — 0, 2; ICrowther 2007|), in 
general the data do not allow for a more accurate determination 
of k, but only to distinguish between these two cases. 

It should be stressed that the intention of our study is not 
to provide ultimate conclusions on the density profiles found for 
individual bursts. Instead, it is meant as a statistical approach us- 
ing bursts with the best available X-ray as well as optical data. 
The questions we want to address are: (1) What is, in a statisti- 
cal sense, the preferred circumburst density profile? (2) What is 
the ratio between events with ISM and with wind profiles? (3) 
What does this tell us about the typical radius of the wind ter- 
mination shock that is expected to exist in the ambient medium 
surrounding a long burst GRB progenitor? 

Throughout the paper we use the convention F v (t) oc t~ a v~P 
for the flux density, where a is the temporal slope and f> is the 
spectral slope. All errors are lcr uncertainties unless noted oth- 
erwise. 



2. Data selection 

2.1. Data gathering 

The optical data were taken from a photometric database main- 
tained and updated by one of the co-authors (D.A.K.). In 
addition, we added data for GRB 090726 from | Simon et all 
(2010), Fatkhu llin etall (120091). lHaislip et all (120091). iKelemenl 
j200%. Landsman & Page! (120091) . ISakamoto et alJ d2009l) and 
IVolnova et al.l d2009l) . The properties of the optical afterglow 
sample, the data gathering and the d e duced SEDs of th e 
afterglows are discussed in iKann et al.l (120061 I2008L l2010h . 
Furthermore, we co mpared their SED results with the work of 
ISchadv et ail (l2007l 12010b . The latter authors added X-ray data 
to model the spectral energy distributions from roughly 1 eV to 



10 keV. Finally, to create a denser light-curve coverage, all non- 
Rc band data in each light curve were shifted to the Rc band. In 
doing this, we used the colours of the SED, assuming no spec- 
tral evolution and omitting data where clear colour evolution was 
evident. 

The X-ray data were retrieved from the Swift data archive 
and the light curves from the Swift light-curve repo sitory (ver- 
sion Februar y 2010) up dated and ma i ntaine d by lEvans etail 
(120071 120091) . Following iNousek et all d2006l) . we reduced the 
data with the software package HeaSoft 6.6.10 together with 
the calibration file version vQ 1 fl Furth e rmore, we app l ied th e 
meth ods detailed inlMoretti et alT(l2005l) . iRomano etafl (|2006), 
and lVaughan et al.l (2006) to reduce pile-up affected data. We ex- 
tracted SEDs at different epochs with approximately 500 back- 
ground subtracted counts, to check for spectral evolution. If 
the properties of the SED (spectral slope and absorption) did 
not evolve, we extracted a new SED from the maximum pos- 
sible time inte rval. In addition, w e used the Tubingen absorp- 
tion model by Wil ms etaD (|2000) and their interstellar-medium 
metal abundance template. The Galactic ab s orptio n was fixed to 
the weighted mean based on Kalb erla et al.l d2005b. We included 
the Ch andra light curve of GRB 05 1221 A from Burrows et al.l 
(200<§. 



2.2. Sample definition and light-curve fitting 

Among all Swift GRBs observed until September 2009, we se- 
lected from our photometric database, which contains long and 
short bursts, those 90 bursts that have an optical and an X-ray 
afterglow as well as a measured spectroscopic or a photomet- 
ric redshift. From these we selected those bursts with the best- 
sampled optical and X-ray afterglow light curves as follows. 

First, we required information on the spectral slope of the 
afterglow in the X-ray band. We derived this information from 
publicly available data. In the optical bands, however, multi- 
band data ar e usually no t available. I n the se cases we u s ed the 
results from IKann et alj (120081 120101) and ISchadv et alJ d2007l 
1201 Oh . with the latter being based on optical-to-X-ray SED fits. 

Second, we required that an afterglow light curve can be 
fitted with a multiply broken power-law and is not dominated 
by flares or bad data sampling. After excluding time intervals 
affected by flares, we fitted the light curves with a smoothly 
broken power-law of the order m (see Appendix [A] for its def- 
inition) with a Simplex and a Levenberg-Marquardt algorithm 
dPress et alj|2007l) . Furthermore, we transformed R c -bmd light 
curves to flux densities with the zero-point definition of Bessell 
dl979l) . For the X-ray regime we followed iGehrels et al.l d2008): 
the flux density, F v x , in //Jy at the frequency v x is then given by 



:4.13x 10 1 



(1 -(3 X )F X 



(lOkeV) 1 ^ -(OJkeV) 1 ^ 



where f3 x is the spectral slope in the X-ray band, F x is the mea- 
sured flux in the 0.3-10 keV range in units of ergcirT 2 s and 
the reference energy E x is given in keV. For all bursts we chose 
the logarithmic mean between 0.3 keV and 10 keV as a refer- 
ence, i.e., £ x =1.73 keV (v x = 4.19 x 10 17 Hz). The numerical 
constant converts the energy in units of keV to a frequency in 
units of Hz and the flux density from erg crrr 2 s _1 Hz 1 to j/Jy. 

Third, once power-law segments had been defined in the 
afterglow light curves, we excluded from further studies those 
bursts where the difference in the late-time decay slopes between 



2 http : //heasarc . gsf c . nasa . gov/ docs/ sof tware/lheasof t 

3 heasarc . gsf c .nasa . gov/docs/heasarc/caldb/swift 
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Table 1. Summary of the considered afterglow models. 



Spectral regime /? x - p opt a x 


- ttopt 


^V,Opt / fy,x 


Spherical expansion 


Vc < V op t < V x 

Vopt < v c < v x 1/2 ± 
Vopt < v x < v c 




1/4 



(v opl /v x ) P 

(Vopt/v x J 


Jet with sideways expansion 


Vc < Vopt < V x 

Vopt < v c < v x 1/2 
Vopt < v x < v c 







(vopt/Vx) '' 

-(p-D/2^,/2 -1/2 
"opt "x "c 

v v p,/v x J 


Jet without sideways expansion 


Vc < Vopt < V x 

Vopt < v c < v x 1/2 ± 

Vopt < V x < V c 




1/4 



(v pt/v x ) 
(V pt/V x j 



Notes. Theoretical differences in the spectral and temporal slopes as 
well as the flux-density ratio depend on the position of the cooling 
frequency, v c , with respect to the optical and the X-ray band (y op t,v x , 
respectively; valid for an electron index , p, of larger than 2; e.g., 
Zhang & Meszaros 2004, Panaitescu 2007). Depending on the circum- 
burst medium the positive (ISM) or negative (wind) solution for a x -a Bpt 
applies. 



the optical and the X-ray band could not be explained within the 
framework of the fireball model (Table [TJ; i.e., the difference in 
decay slopes, a x - a opt , was larger than 1/4 within 3cr. Because 
of this, a strict criterion for the beginning of the late-time evolu- 
tion of an afterglow cannot be given. Evidence for an observed 
canonical light curve does not exist in th e optical dKann et al.l 
2008, l2010t iPanaitescu & Vestrandll2010l). but may exist in the 
X-ray s (e.g. jNousek et al.ll2006t IZhang et al.ll2006t lEvans et all 
2009). Therefore, as an operational definition an afterglow is in 
its late-time phase when its temporal and spectral evolution can 
be explained by the fireball model from a particular time after 
the corresponding GRB. 

After applying these selection criteria, about half of the 90 
bursts had to be rejected owing to bad sampling, poor data qual- 
ity, flares, and other peculiarities. Another quarter had to be re- 
jected because of \a x - a opt \ > 1 /4 within 3<x. All data were then 
corrected for host extinction in the optical, Ay OSt , if needed, and 
for Galactic and host absorption in the X-ray band, A^ ost . 

In total, 27 afterglows passed these selection criteria. The 
sample consists of 25 long and one short GRB (GRB 05 1221 A), 
and one controversial event (GRB 060614) in terms of the 
short/long classification scheme. Their input data are sum- 
marised in Table lB.ri and lB.2l 

3. Results and discussion 

3.1. The circumburst medium density profile 

3.1 .1 . Identifying the circumburst medium 

As a first step in the identification of the circumburst medium, 
we defined nine possible spectral and dynamical regimes (Table 
|2J. The spectral regimes separate according to the position of the 
cooling frequency with respect to the observer frame, while the 



Table 2. The closure relations combining the temporal 
decay slope a and th e spectral slope & (adopted from 
Zhang & Meszaros 2004 and Panaitescu 2007; valid for p > 2). 



Spectral Closure relation a (J3) 

regime ISM Wind 

Spherical expansion 

v c > v 3/3/2 [Sla] (3/3 + l)/2 [Sib] 

v c <v (3/?-l)/2[S2] 

Jet with sideways expansion 

v c >v 2/J+l[Jl] 
v c < v 2/3 [12] 



Jet without sideways expansion 

v c > v (6/3 + 3)/4 [j la] (3/3 + 2)/2 [j lb] 
v c < v (6/3 + l)/4 [j2a] 03 + 5)/4 [j2b] 



Notes. Abbre viation used here for a certain model is given in brackets 
(adopted from Panaitescu 2007; (S, J, j)=(spherical expansion, jet with 
sideways expansion, jet without lateral spreading), (1, 2)=(v < v c , v > 
v c ), (a, b)=(ISM, wind)). Entries extending over two columns are valid 
for both media. 

dynamical regimes distinguish between a spherical and a jetted 
evolution. 

To derive the most probable density profile into which a GRB 
jet propagates, we proceeded in the following way. The main 
criterion was that the result agrees with the optical as well as 
with the X-ray data. In doing so, we first analysed the closure 
relations for the nine models (Table |2]i in the optical and X-ray 
bands and selected only those relations (models) which were ful- 
filled within 3cr. Second, we computed the difference in the de- 
cay slopes, a x - a opt , and distinguished between the models ac- 
cording to the three possible cases -1/4, 0, + l/4(Table[TJl, again 
within 3<x. Third, if possible, we took into account the difference 
in the spectral slope, /3 X -/3 opt , which is either or 1/2 (see Table 
[TJ. If this criterion could be applied, we required that it is ful- 
filled within 3<r. In addition we required that the Icr uncertainty 
in the spectral or temporal decay slopes was less than 0.2. 

The final results of the light curve fits of the 27 bursts con- 
sidered are presented in Fig. IC.ll In addition, we plot here 
the observed flux-density ratio (f v.opt/^v.x) (0 (middle panels in 
Fig. IC.lt as a function of time. According to Table [TJ the ex- 
pected flux-density ratio is only allowed to take a certain value 
depending on the spectral and dynamical regime. We checked if 
the observed flux-density ratio, \F v ,opt/F VyX ^ (t) (middle panels in 
Fig. IC.ll ). agreed with the model(s) that successfully passed the 
previous three criteria. The allowed parameter space of the flux- 
density ratio of all considered models (Table[TJ) is shown as grey 
box in the middle panels in Fig. IC.ll The upper and lower bound- 
ary always refer to v c < v opt with F Vi0pt /F v , x = (v opt /v x y p/2 and 
to v c > v x with F Vt0pt /F VtX = (v opt /v x )~ ( - p _1) / 2 . This criterion came 
into play when we were unable to distinguish whether the cool- 
ing frequency was redward of the optical band or blueward of the 
X-ray band, in other words when the flux-density ratio agreed ei- 
ther with the lower or upper boundary in Fig. IC.ll 

3.1.2. Ensemble properties 

Combining these different criteria, we found that for about 60% 
(16/27) of all studied cases the cooling break was between the 
optical and the X-ray band (Table lB~3l . i.e., the difference in the 
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decay slope is either + 1/4 or -1/4 (for a constant density pro- 
file and a free wind medium, respectively). In total, we could 
identify the circumburst density profile (ISM or wind) for 25 of 
the 27 investigated afterglows (Table iBTJb . However, we identi- 
fied a wind medium for only six events (GRBs, 050603, 07041 1, 
080319B, 080514B, 080916C and 090323). The other 19 bursts 
were consistent with an ISM profile except for GRBs 05 1221 A 
and 060904B. 

Our procedure to find physical descriptions of light-curve 
segments, in other words identifying the spectral regime and the 
circumburst medium, allowed us to find descriptions of more op- 
tical and X-ray afterglows than are presented in the literature 
(Table IB.3b . Our results usually agree with the literature (for 
references see Table IB. 3k they only differ for GRBs 070802, 
080721, 090323, and 090328. The contradiction in the latter 
three bursts is due to the size of the optical data se t. We used the 
maxim um publicly availabl e data set in contrast to lStarling et al.l 
(2009) (GRB 080721) and lCenko et ail (1201 Oh (GRBs 090323 
and 090328). 

Four bursts are of p articular interest in our sample: (i) GRB 
051221A (z - 0.54 6; ISoderberg etail 120061) is a short burst 
dBurrows et al.l I2006T) . i.e., m ost likely it originated from the 
merge r of two co mpact objects (jBlinnikov et al.ll98"4tlPacz vhski 
1986; iGoodmanl fl986t lEichler et al.l ll989l). Unfortunately, we 
cannot discern between a wind or an ISM profile because v c < 
v x . The cooling frequency was below the optical bands so that 
neither the optical nor the X-ray data can be used to reveal the 
circumburst medium except during a post-jet break phase with- 
out lateral spreading. 

(ii) GRB 060614 (z = 0.125; iDella Valle et"ail l2006bh 
is a quite controversial event that does n ot easily fit into 



the classical short/long c lassificat ion schem e (Delia Valle et al. 
2006 a]: iFvnbo et all 120061; [Gal- Yam et all l2006t iGehrels et al l 
2006: lMangano et al.ll2007l: Izhang et al.ll2007l:lKann et al.ll2 008 ; 
Zhang et all 120091) . The difference in the spectral slope, /3 X - 



p opt = 0.00 + 0.11 (Table |B3]l, rules out that the cooling fre- 
quency lies between the optical and X-ray bands. Furthermore, 
we did not find evidence for a wind medium. The optical and 
X-ray data exclude a wind medium with high confidence. The 
deviation between the observed and predicted temporal decay 
slope is < 1.3cr and > 3.5<x for an ISM and a wind medium, 
respectively. In our sample this burst belongs to a small number 
of cases where the flux density ratio could be used to distinguish 
between v c > (v op t, v x ) and v c < (v opt , v x ). The small observed 
flux density ratio agrees with v c > (v opt , v x ). 

(iii) GRB 060904B has a well-defined light curve (Fig.lCTI 
Table IB. 2\ and SED (Table iBTTT i in the optical and X-ray bands, 
respectively. The difference in the decay slopes, a x - a opt = 
0.20 + 0.04 (Table |S7J), favours an ISM profile with the cool- 
ing frequency lying between the optical and the X-ray bands. 
The difference in the spectral slopes, however, does not support 
this scenario, /3 X - /3 opt = 0.00 + 0.16 (Table E3J. Rather, both 
afterglow components seem to be in the same spectral regime. 
Therefore, we could not find a consistent description of the op- 
tical and X-ray afterglow. The main reason could be t hat the op- 
tical d ata are mainly based on preliminary data (see iKann et all 
120101 for details on the data gather ing). 

(iv) GRB 0803 19B (z = 0.937: lD'Elia et al.l2009l) is the only 
burst in our sample with a photometrically detec ted supernova 
component (Bloo m et al 1120091; iTanvir et alj|2008l) . An ISM pro- 
file can be ruled out with very high confidence. The difference in 
the spectral slopes, /3 x -/3 opt = 0.48+0.12 (Table lB~3l . favours the 
cooling break to be in between the optical and the X-ray bands. 
The deviation between Aa b s - a x - a opt and Aa pre di cte d is 1 <x 



for a wind medium but 9<x for an ISM profile. The closure re- 
lations support this finding. The deviation between the observed 
and predicted optical decay slope is . 1 cr for a wind m edium and 
4.6cr for an ISM profile (see also lRacusin et al.ll2008h . 

3.1.3. The electron index 

Finally, the identification of the light-curve segments allowed 
us to derive the electron index, p. It is shown in the top pan 
els in Fig. IC. 1 J for every burst. In agreement with other stud- 
ies (e . g.,[P anaitescu & Kumar 2001; Shen et al. 2006 : IZeh et al 
2006; Starling et al.l 120081: ICurran et all l2009t iGhisellini eta! 



2009: ICurran et al.ll2010l) we find that the distribution extends 
from p ~ 2 to p ~ 3, indicating that p is no universal value. 

3.2. Constraining the wind - ISM transition zone 

Once we had identified the circumburst density profiles, we 
could investigate if there is observational evidence for the po- 
sition of the wind termination-shock radius, where the density 
profile changes from that of a free wind (k = 2) to a density pro- 
file with k - (independent of whether it is the shocked wind or 
the ISM; see lPe'er & Wiiersl2006l:lvan Marie et al]2007h . While 
the theory of a blastw ave crossing such a d ensity discontinuity 
has been worked out dPe'er & Wiiersl 120061) . finding the corre- 
sponding observational signature is difficult. Even though data 
on several hundred afterglows exist, they do not provide this 
kind of informatio n in a convincing way (IStarling et al.l 12008: 
ICurran et al. 1 120091 This leaves open the question of observa- 
tionally determining the wind termination-shock radius. 

Here we cannot determine the radius of the wind termination 
shock for any (long) GRB progenitor either. However, we can 
characterise its position in a statistical sense. The light curves 
in Fig. IC.ll show F opt /F x as a function of time. Its first log- 
arithmic derivative, (a x - a opt ) (f), reveals either the maximum 
time up to which a wind profile is identified or it reveals the 
minimum time after which an ISM profile agrees with the data. 
The function (a x - a opt ) (f) is a smooth function in time in con- 
trast to the fit values, because we used a smoothly broken power 
law to describe the light-curve evolution. Table[3]summarises the 
time intervals for which the asymptotic values for a x - a opt were 
reached, i.e. -1/4, 0, or 1/4 within l<x. 

In the observer frame, t he radius of the fireba ll propagating 
into a free wind medium is dChevalier & Li 2003) 



R(t) = 1.1 x 10 1 



2f£ is 



iso.52 



(l+z)A* 



1/2 



cm , 



(1) 



where t is measured in units of days, Ei so is the isotropic equiv- 
alent energy in units of 10 52 erg and A* is defined via A = 
M w /47tv w = 5 x 10 11 A*gcm , with M w being the mass-loss 
rate, and v w the wind velocity. The quantity A* refers to a mass- 
loss rate of M w = 10~ 5 Mq yr~' and a velocity of the stellar wind 
of v w = 10 8 cms -1 . 

Using the aforementioned approach and fixing for simplic- 
ity A* = 1, Table [3] summarises the deduced upper and lower 
limits on the wind termination-shock radius. The distribution 
extends over three orders of magnitude (from a 10~ 3 pc to 1 
pc). This spread is partly due to the strong correlation with Ei S0 

1 /2 

(R(t) oc EJ ) and is thus related to the spread in the energy re- 
leased during the prompt emission in gamma-rays (width a 4.3 

— 1/2 

dex). On the other hand, since R(t) scales with A . , A* would 
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Table 3. Constraints on the termination-shock radii. 



GRB 



'start 

(ks) 



*end 

(ks) 



log E isi 
[erg] 



R T 
(pc) 



Constant density medium (n(r) oc r°) 



050801 


0.31 


105.70 


050820A 


42.71 


411.00 


05 11 09 A 


25.87 


264.60 


060418 


4.40 


535.80 


060512 


10.50 


287.50 


060614 


51.59 


1497.00 


060714 


78.49 


286.10 


060908 


0.92 


82.14 


070802 


10.90 


88.78 


07081 OA 


1.95 


34.41 


080710 


23.10 


348.50 


080721 


29.46 


1258.00 


08 1203 A 


12.07 


301.00 


090102 


77.57 


262.80 


090328 


57.89 


923.30 


090726 


3.61 


10.00 


090902B 


45.64 


1168.00 


090926A 


325.80 


1798.00 



51 51 +034 

-0 12 

53.99!:' 
52.87+ : 

53 15 +8 -" 

50 30+"- 40 

51 40 +005 

-Q04 

52 89 + "- 3 ° 

52.82 ± 0.20 

51 70 +031 
-QQ9 

51 96 +m 

-815 

51.90+"' 

54 09+ 03 

-Q04 

53 54 + ° A * 

JJ ' J -0.15 

53.30 ± 0.04 
52.99 ± 0.01 

52 26 +033 
54.49 ± 0.01 
54.27 ± 0.02 



< 0.001 « 

< 0.18*$ 

< 04+° 01 

< 027 +8 o ° 4 

< o 002 r° 0012 

<- U.UUil _o.oOQ2 

< 0.018 ±0.001 

< 069+ 0029 

< 0.008 ± 0.002 

< 007 +0003 

-o.no l 

< 004 1 + ° 0002 

u ' wtl -0.0007 

< 0.02 ± 0.01 

< 0.17 ±0.01 

< 0.06 ± 0.01 

< 0.13 ±0.01 

< 0.098 ± 0.001 

< 007 +0 003 
<- u.uu/ _ 0001 

< 0.383 ± 0.004 

< 0.76 ± 0.02 



Free stellar wind (n(r) oc r 2 ) 



050603 
070411 
0803 19B 
0805 14B 
0809 16C 
090323 



39.21 
92.09 
52.76 
37.30 
97.49 



196.50 
522.70 
555.70 
174.70 
377.70 



53.79 ± 0.01 > 0.305 ± 0.004 



53 00+ 26 
54.16 ±0.01 

53 42+ 03 

JJ -^ -0.04 
za 4 q+0.08 
J1, -0.10 



> 20+ 07 

> 1.11 ±0.01 

> 0.22 ± 0.01 

> 0.80^ 



73.05 1079.00 54.61 ±0.01 > 1.68 ±0.02 



Unclear cases 


050922C 

051221A 25.92 
060904B 3.58 


228.80 
162.80 


52 98 +oos 

-QQ5 
51 41+°- ifc 

J - ^ -0.33 

51 71+ 013 

J1 - ' -0.21 



Notes. Time intervals deduced from the lower panels in Fig. IC.ll if 
the first logarithmic derivative of F opt /F x points to either an ISM pro- 
file or a free wind, i.e., once the asymptotic value +1/4 or -1/4 was 
reached (within t he error bars). With the exception for G RB 0805 14B 
dRossi et aUfcOOSft and GRB 090726 dButler et alj|2010h. the iso tropic 
equivalent energies, E lso , were taken from Kann et al. (2008, 2010). The 
forth column was calculated based on Eq. [T] assuming A, = 1, with 
t = faun (Fig. IC.lb for ISM profiles and t = F em3 (Fig. IC.lt for wind 
profiles. 



have to vary by a factor of 100 in the right way to reduce the 
width of this distribution by only a factor of 10. 

Is the width of the distribution we have found for the up- 
per limits on the wind termination shock rea sonable? Based on 
numerical wind models of Wolf-Rayet stars, iFrver et all (2006) 
found that for A* = 1 the radius of the termination shock is ap- 
proximately given by 



fiT = (»iSM/100ctrr 3 r 1/2 pc. 
ISM densities of the order of 10 6 ™""" 3 



(2) 



cm J are then required to re- 
duce Rj to 0.01 pc. These large-scale gas densities are rather 
unique and only typical for dense cores of mol ecular cloud s 
(with the Rho Ophiuchi Cloud as an example, e.g., Klose 1986). 

Observationally, gas densities could only be derived for a few 
bursts because of the lack of radio data. So far, the highes t values 
are about 600 cm 3 dFrail et al.ll2006t iThone et al.ll2010h. while 
the n ominal value is about a factor of 100 smaher dFrail et all 
2006). These measurements do not necessarily rule out the 




10" icy 1 
Radius R(t) (pc) 

Fig. 1. Shown here are lower (to-the-right pointing triangles) and 
upper (to-the-left pointing triangles) limits on the position of the 
wind termination shock based on Eq. [TJ assuming A* = 1 in all 
cases (TableQ. Note that GRB 060614 is a much debated burst 
(see Sect. 13.1. Ti l. The step curves are the cumulative distributions 
to the lower limits up to which a wind profile is identified and 
upper limits after which a constant density medium (ISM) agrees 
with the data. 



model of IFrver et aD (120061) . Successful radio observations 
might have picked out a certain class of GRBs. Furthermore, 
radio observations are very challenging at early times because 
the brightness of the afterglow is increasing in the radio bands 
while t he afterglow is already de caying in the optical and X-ray 
bands dZhang & Meszarosll2004l) . Thus, it is difficult to extract 
information on the direct vicinity of the progenitor. On the other 
hand, if the average particle density of the circumburst medium 
is about 1-10 crrT 3 , then additional mechanisms are required to 
bring the wind term ination-shock radius closer to the star (e.g., 
Ivan Marie et al]|2006l) . 

The other possibility to reduce Rj is to decrease A*, since 
R T oc A\l 2 dChevalier et al.ll2004l) . For example, for A* = 0.01 
all data points in Fig. [TJ would shift along the /^-coordinate to 
higher values by a factor of 10 (Eq.Q}, while Rj would decrease 
by the same factor. In this case, lower circumburst gas densities 
would be required. This touches upon the question on how small 
A* can be. Studies of W R stars do not f avour values of less than 
0.01 in polar directions (Eldridge 2007, his Table 1). Moreover, 
observations of ne arby WR stars do not sh ow evidence for these 
low values either (iNugis & LamersluOOOl) . However, the major- 
ity of nearby WR stars are surely not seen pole-on, in contrast to 
GRBs. Therefore, it is difficult to decide if these observational 
constraints on WR stars can be applied to GRB progenitors. 

On the other hand, the six GRBs that favour a free wind 
medium (Table IB.3b have large lower limits on the wind 
terminatio n-shock radius ( Table |3). This matches theoretical 
models by Ivan Marie et al.l d2007lEo08l) . which allow R T to ex- 
tend up to several parsecs. 

The separation between the lower and upper limits for wind 
and ISM-profiles, respectively, on the wind termination-shock 
radius could be even larger. Refining the lower and upper limits 
is difficult, however. The lower limits, which are deduced from 
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fend (Fig. IC.ll Tabled, depend on the observing strategies due to 
the brightness of the afterglow and the brightness of the underly- 
ing host galaxy in the optical bands. On the other hand, the upper 
limits, which are deduced from f start (Fig. IC.ll Table |3), can be 
affected by additional radiation components at early times. 

4. Summary and conclusion 

After applying several selection criteria (closure relations, the 
differences in the spectral and temporal slopes, and the flux den- 
sity ratio, F opt /F x ) we selected the best-sampled Swift GRBs 
with well-observed optical as well as X-ray afterglow data from 
June 2005 to September 2009. Altogether 27 bursts entered our 
sample, which was used to investigate the density profile of the 
circumburst medium (constant density medium or free wind), 
including one short burst (GRB 05 1221 A) and one controversial 
event in terms of classification (GRB 060614), which success- 
fully passed our selection criteria among all bursts. The other 25 
events are classified as long bursts without doubt. 

Combining optical with X-ray data is advantageous because 
optical data usually allow for a more precise determination of 
the temporal decay slope of an afterglow, while X-ray data can 
in general be used to extract the SED. Combining both emission 
components substantially improves our capability to distinguish 
between an ISM and a wind medium. Thereby, we concentrated 
on the late-time evolution, i.e., times when the proper afterglow 
is not affected anymore by flares and additional radiation com- 
ponents (e.g., the reverse shock, central engine activity). 

Our study shows that only six of the 25 long bursts (24%) 
investigated here (GRBs 050603, 070411, 080319B, 080514B, 
0809 16C and 090323) showed evidence for a free wind medium 
at late times. In the other cases (76%), except for the short burst 
GRB 05 1221 A and 060904B, the blastwaves were propagating 
into a constant density-medium. In particular, the controversial 
burst GRB 060614 favours an ISM profile. This is not in dis- 
agreement with a massive-star origin as our result for long bursts 
indicates. 

In addition, we were able to set limits on the wind 
termination-shock radii of the corresponding GRB progenitors. 
Only 24 of 27 bursts (Table [3} had good enough data to per- 
form this analysis. Fixing the relative mass-loss rate to A* = 1, 
the distribution we deduced covers three orders of magnitude. 
We find a tentative grouping into (long) GRB progenitors with 
comparably small and comparably large termination-shock radii. 
Whether this points to two distinct populations of (long) GRB 
progenitors or if this is a selection effect, remains an open issue. 
At least theoretically it is well possible that the long burst pop- 
ulation splits int o single star proge nitors and those belonging to 
a binary system dGeorgy et al.ll2009l) . Further observational data 
are required to reveal a potential binary nature of the long burst 
progenitors. 
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Appendix A: Smoothly broken power law of the 
order m 

The equation for a smoothly broken power-law of the order m, 
F v , m (t), was derived by recursion in the following way. Let us 
assume the function F vm (f) consists of m power-law segments 
connected by (m — 1) breaks. To add an additional power-law 
segment F Vy „ 1+ \(f), we first normalised the new power-law seg- 
ment to the previous one, F vm {t), at the break time tb, m - 



F v ,m+i(t) - F Vim (tb jn ) 

\tb,m 



(A.l) 



Here a m +\ is the slope o f segment (m +1). Second, we followed 
Beuermann et al ] (119991) and introduced a smoothness parameter 
n m so that the smoothly broken power law of order the (m + 1 ) 
takes the form 



F Vtm + 1 (t) = (F-%(t) + F- v % +1 (tj) 



(A.2) 



If the light-curve consists of m segments, both steps (adding and 
smoothing) have to be performed (m - l)-times. 

For exa mple, let us derive the e quation for a smoothly broken 
power-law dBeuermann et al.lll999l) . In this case m = 2, thus the 
function consists of two power-law segments connected by one 
break at the time tb,\. The initial function is a simple power law 
F v ,\(t) - Ct~ ai . First, the second power-law segment, F v ^, has 
to be connected to the first one at the time t/,,i (step lA.lt 



F v ,2 - F % \ (ffo.i) — 
\tb,i 



Ct 



h.\ 



tb,\ 



Second, the transition has to be smoothed by weighting both 
functions at the point of intersection (step lA.2t 



F v , 2 = (F-£(t) + F- v £(t)) 



-l/Hl 



t 



■l/m 



ai n\ \ — 1/n 



This leads to the equation found by Beuerm ann et alJ (1 19991) for 
a smoothly broken power-law. Repeating both steps leads to a 
smoothly broke n power-law of the order 3 (double smoothly bro- 
ken power-law: iLiang et al.l l2008). Thus, looping (m - l)-times 
over both steps results in a smoothly broken power law of the 
order m. 



Appendix B: Tables 
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Table B.l. Properties of the afterglow SEDs in the optical and X-ray bands of the 27 bursts that entered our sample. 



GRB 


z 


Popt 


Dust 
Model 


A host 
V 


ft, 


i mid 

(ks) 


References 


050603 


2.818 










0.96 ±0.12 


X: 44 


z: 1; Opt: 2 


050801 


1.560 


V.Oy ± 


U. J^- 


SMC 


0.30 + 0.18 


0.89 ± 0.08 


Opt: 86; X: 79 


z: 3; Opt: 2 


050820A 


2.615 


O 77 -i- 
U. 1 L ± 


U.Uj 


SMC 


0.07 ± 0.01 


1.09 ±0.03 


Opt: 86; X: 37 


z: 4; Opt: 2 


050922Ct 


2.199 


1 07 




MW 


0.17 ±0.05 


1.07 ±0.05 


Opt/X: 20 


z: 6; Opt/X: 5 


05 11 09 At 


2.346 






SMC 


< 0.10 


0.90 ± 0.04 


Opt/X: 5 


z: 7; Opt/X: 5 


051221A 


0.546 










0.95 ±0.11 


X: 3 


z: 8; Opt: 9 


060418 


1.490 


O f^Q -t- 


n 1 1 


LMC 


0.20 ± 0.08 


0.98 ±0.15 


Opt: 86; X: 7 


z: 10; Opt: 2 


060512 


0.443 










1.02 ±0.11 


X: 23 


z: 11 


060614 


0.125 


U.o 1 ± 


U.Uo 


SMC 


0.05 ± 0.02 


0.81 ±0.11 


Opt: 67; X: 36 


z: 12; Opt: 13 


0607 14f 


2.711 


n Q7 




SMC 


0.46 ±0.17 


0.92 ±0.10 


Opt/X: 5 


z: 14; Opt/X: 5 


060904B 


0.703 


1 1 1 -L. 

1. 1 1 ± 


n 1 n 


SMC 


0.08 ± 0.08 


1.11 ±0.12 


Opt: 86; X: 2 


z: 4; Opt: 2 


060908 


1.884 


U.JU ± 


U.Uj 


SMC 


0.00 


0.94 ± 0.08 


Opt: 86; X: 2 


z: 4; Opt: 2 


070411 


2.954 










1.19 ±0.14 


X: 8 


z: 4 


070802t 


2.454 


0.61 




MW 


1.20 ±0.12 


1.11 ±0.05 


Opt/X: 2 


z: 4; Opt/X: 15 


0708 10A 


2.170 










1.14 ±0.16 


X: 2 


z: 16 


0803 19B 


0.937 


0.50 ± 


0.07 


SMC 


0.15 


0.98 ±0.10 


Opt/X: 170 


z: 17; Opt/X: 18 


0805 14Bf 


1.800 


0.63 ± 


0.02 


SMC 


0.00 


1.13 ±0.13 


Opt: 43; X: 33 


z/Opt/X: 19 


0807 10t 


0.845 


1.01 




SMC 


0.00 ± 0.00 


1.01 ±0.01 


Opt/X: 27 


z: 4; Opt/X: 20 


080721 


2.591 










0.99 ± 0.04 


X: 4 


z:4 


0809 16Ct 


4.350 


0.49 




SMC 


0.00 


0.49 ± 0.34 


Opt: 79; X: 67 


z/Opt/X: 21 


08 1203 A 


2.050 










1.06 ±0.07 


X: 22 


z: 22 


090102 


1.547 


0.74 ± 


0.22 


SMC 


0.12 ±0.11 


0.77 ± 0.04 


Opt: 86; X: 8 


z: 23; Opt: 2 


090323 


3.568 


0.65 ±0.13 


SMC 


0.14 ±0.04 


0.95 ±0.13 


Opt: 97 ks, X: 92 ks 


z: 24; Opt: 25 


090328 


0.735 


1.17 ± 


0.17 


SMC 


0.18 ±0.13 


1.10 ± 0.16 


Opt: 86; X: 130 


z: 25; Opt: 2 


090726 


2.710 










1.45 ±0.15 


X: 17 


z: 26 


090902B 


1.822 


0.73 ± 


0.13 


SMC 


0.05 ± 0.05 


1.02 ±0.11 


Opt: 86; X: 98 


z: 24; Opt: 2 


090926At 


2.106 


1.04 




MW 


< 0.10 


1.04 ±0.08 


Opt/X: 250 


z:24; Opt/X: 27 



Notes. For every GRB we list the redshift, z, the properties of the optical SED (spectral slope {3 opt , host extinction in the V-band rest frame, Ay OSt , 
and the corresponding extinction law), the spectral slope in the X-rays, /3 X , and the mean time after the burst, f raid , at which the SED was extracted. 
The events for which we used joint optical-to-X-ray SEDs are marked with a 'f. In these cases the difference in the spectral slope, /3 X - [} opi was 
either fixed to 1 /2 or 0. Except for GRB 0805 14B, these estimates have only one error estimate for /3 opt and f3 x together due to the fitting procedure 
described in the associated papers. 

Refer ences, m iBerger & Becked J2005h: (2) iKann et all d201Ch: (3) Ide Pasauale efaD j2007h: (4) iFynbo et alj d2009h: (5) ISchadv et alJ d2010h : 
(6) iPiranomonte et alj d2008h: (71 lOuimby et alj j2005h : (8) ISoderberg et alj j2006l) ~( 9') IKann et alj d2q08r);~ (10) iProchaska et alJ j200d) : (11) 
Bloo m et alj d2006h : (12)" lDella Valle et alJ j2006bl); 7l3) iMangano et al.N2007h: (14) Uakobsson et al.l j2006h: (15)TKruhler et alJ 12008): (16) 
Thon e et al.U2007l): (17) i D'Elja et al.l d2009T): (18) iRacusin et all d2008l) ; (19) iRossi et alj j2009l) ; (20) iKriihler et alj j2009h ; ( 21) breiner etaL 
( |2009if ): (22 ) iKuin et al l 00091) : (23) fde Ugarte Postigo et alJd2009h : (24) ICenko et alJd2010h : (25) iMcBreen et alJd2010T) : (26) iFatkhullin e~ 
(2009): (27) feau et alj d201Ch 
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Table B.2. Light-curve parameters of the late-time optical and X-ray afterglows of the 27 bursts that entered our sample. 



GRB 



Band 



Light-curve parameters 
(ks) «(ks) ' or, a 2 



Overlapping time interval 

fbjel (ks) f 6taft (ks) fend (ks) 



050603 

050801 

050820A 

050922C 

05 11 09 A 

05 1221 A 

060418 

060512 

060614 

060714 

060904B 

060908 

070411 

070802 

07081 OA 

0803 19B 






33.33 


196.67 




1.97 ±0.06 


X 


39.21 


577.04 




1.67 ±0.05 





0.25 


106.27 


1.19 ± 0.01 




X 


0.22 


448.10 


1.29 ± 0.08 







23.89 


612.10 


1.04 ± 0.01 




X 


4.64 


3676.50 


1.19 ± 0.02 


1.73 ±0.16 





4.13 


606.01 


1.47 ± 0.04 




X 


2.02 


95.92 


1.37 ±0.03 







17.52 


265.20 


1.03 ± 0.06 




X 


3.46 


1240.80 


1.20 ± 0.01 







4.68 


445.12 


0.96 ± 0.03 




X 


15.17 


2271.92 


1.03 ± 0.01 


1.99 ±0.17 





0.12 


871.05 


1.192 ±0.002 




X 


4.40 


537.08 


1.52 ±0.05 







7.02 


715.56 


0.80 ± 0.03 




X 


0.24 


288.61 


1.12 ±0.05 







51.59 


1700.19 


1.05 ±0.04 


2.45 ± 0.05 


X 


41.23 


1500.57 


1.07 ±0.11 


2.31 ±0.11 





47.07 


286.77 


1.42 ± 0.18 




X 


3.07 


1029.03 


1.23 ±0.03 







3.52 


163.13 


1.16 ± 0.02 




X 


3.58 


391.01 


1.36 ± 0.03 







0.14 


82.42 


1.03 ± 0.01 




X 


0.71 


483.21 


1.49 ±0.07 







92.09 


523.12 


1.29 ± 0.06 




X 


0.50 


664.68 


1.12 ±0.02 







10.90 


88.88 


0.90 ±0.16 




X 


6.68 


316.52 


1.17 ±0.09 







1.95 


103.36 


1.30 ±0.10 




X 


1.55 


34.46 


1.29 ±0.07 







0.8 


1060 


1.237 ±0.002 




X 


42.66 


2559.17 


1.04 ±0.05 


2.68 ± 0.36 



< 39.2 

> 106.3 

640.2 ± 138.7 

> 95.9 

> 265.2 

318.5 ±38.6 

> 537.1 

> 288.6 

113.1 ±2.7 
120.0 ± 12.1 

> 286.8 

> 163.1 

> 82.4 

> 523.1 

> 88.9 

> 34.5 

953 ± 137 



39.2 
0.25 
23.9 

4.1 
17.5 
15.2 

4.4 

7.0 
51.6 
47.1 

3.6 
0.71 
92.1 
10.9 

1.9 
42.7 



196.7 
106.3 
612.1 
95.9 
265.2 
445.1 
537.1 
288.6 
1500.6 
286.8 
163.1 
82.4 
523.1 
88.9 
34.5 
1060 



Notes. Columns 3 to 7 give the time interval, f'^ t , 4? p, the pre- and post-jet break decay slopes, 0-1,0-2, and the observed jet break time, f^e,. The 
optical and X-ray light curves were independently fitted for every GRB (Sect. l2.2V Owing this, two break times are given for GRB 060614. Within 
the errors the break was achromatic, making it a good candidate for a jet break. For the other GRBs we state the upper or lower limits on the jet 
break time with respect to the identification of the dynamical regime shown in Table IB. 31 The overlapping time interval of the late-time optical 
and X-ray afterglow is shown in the last two columns, f The optical afterglow light curve of G RB 080721 shows an additional shallow break at 
(129 ± 84) ks. The difference in the pre- and post-break decay slope is sa 0.25 in agreement with Kan n et al] d2010h . which is typical for a cooling 
break. This break is not a jet break. 
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GRB 


Band 




Light-curve parameters 




Overlapping time interval 


fstarl (ks) 


fend (KS) 




<1'2 


fbjet (ks) 


fstait (ks) 


fend (ks) 


0805 14B 


O 

X 


37.15 
37.30 


174.79 
217.27 


1.64 + 0.06 
1.54 ±0.14 




> 174.8 


3 1 .i 


1-7/1 O 

1 /4.S 


080710 


O 

X 


9.97 
11.34 


353.11 
349.31 


1.57 ±0.02 
1.56 ±0.09 




> 349.3 


i 1.3 




08072 If 


o 

X 


0.17 
29.46 


2641.08 
1256.83 


1.22 ±0.01 
1.50 ±0.03 


1.46 ±0.08 


> 1256.8 


29.3 


1236.8 


0809 16C 


o 

X 


97.49 
65.03 


377.94 
1306.33 


1.40 ±0.10 
1.29 ±0.08 




> 377.9 


9/. 5 


3 / l.y 


081203A 


o 


6.60 


301.54 




1.72 ±0.01 




o n 


3U1.3 


X 


8.91 


345.00 


1.13 ±0.01 


1.93 ±0.06 


8.9 ± 0.9 


8.9 


090102 


o 

X 


15.76 
0.95 


263.68 
688.48 


1.50 ±0.03 
1.43 ±0.02 




> 263.6 


15.8 


263. / 


090323 


o 

X 


96.81 
70.47 


1142.30 
1084.77 




1.88 ±0.01 
1.56 ±0.10 


< 96.8 


yb.a 


i no a o 
11)54.8 


090328 


o 

X 


57.28 
0.15 


1007.14 
924.60 




1.78 ±0.04 
1.68 ± 0.09 


< 57.3 


57.3 


924.6 


090726 


o 

X 


1.40 
3.61 


10.00 
66.43 


0.97 ±0.10 
1.34 ±0.04 




> 10.0 


3.6 


10.0 


090902B 


o 

X 


45.64 
45.21 


1171.46 
1456.69 


0.97 ± 0.02 
1.33 ±0.03 




> 1171.5 


45.6 


1171.5 


090926A 


o 

X 


254.24 
51.49 


2070.84 
1803.53 


1.74 ±0.02 
1.53 ±0.07 




> 1803.5 


254.2 


1803.5 



Table lB.2l — continued 
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Table B.3. Identification of the light-curve segments and the circumburst medium. 



GRB 


Closure relation 




A- 


Popt 


Conclusion 


Literature 
ID Reference 


UJUUUJ 


O: — 

X: jl/2a,j2b 


— o ^o + o or 

— U.JU IE u.uo 






U. JlD 

A. JZD 


Wind 






UJUOU 1 


O: 

X: Sla 


fl 10 + f) 08 

U.lu I \J.\JO 




+ 35 


0:Sla 

Y- S 1 fl 


ISM 


O: Sla 

Y- S 1 n 


1 


050820A 


O: Sla 
Xi: S2 
X 2 : J2,j2a/b 


0.15+0.02 


0.37 


±0.04 


O: Sla 

v . CO 
A] . oZ 

X 2 : j2a 


lolVl 


O: Sla 

A[ . oZ 

X 2 : — 




L 




O: Sla 
X: Sla 


-0 10 + OS 


00 


+ OS 


X: Sla 


ISM 






0S1 10QA 


O: 

X: Sla, S2 


1 7 + Ofi 

U. 1 / IE VJ.vJU 


so 


+ 04 


O: Sla 
X: S2 


ISM 








O: — 


0.07 ± 0.03 














051221A 


Xi:Sla, S2 

X 2 : J2,jl/2a,jl/2b 






X,: S2 


— 


X x : S2 


3 


060418 


O: Sla 

X: Sla/b, S2 


0.33+0.05 


0.29 


±0.19 


0:Sla 

V- QO 


ISM 






060512 


O: — 

X: Sla, S2 


0.32 ± 0.06 




— 


O: Sla 

A. oZ 


ISM 






UUUU14 


Oi:Sla, S2 
Xi: Sla, S2 
2 : Jl, jib 
X 2 : Jl/2, jla/b 


0.02 + 0.12 
-0.14 + 0.12 


0.00 
0.00 


±0.11 
±0.11 


Oi: Sla 
X,: Sla 
2 : Jl 
X 2 : Jl 


ISM 


X,: Sla 
X 2 : Jl 


4 


060714 


O: Sla/b, a 2 
X: Sla, S2 


-0.19 + 0.18 


0.00 


±0.10 


OrSla 
X: Sla 


ISM 






060904B 


O: S2 

X: Sla, S2 


0.20 ± 0.04 


0.00 


±0.16 


O: S2 
X: S2 








060908 


O: — 
X: Sla 


0.46 ± 0.07 


0.64 


±0.09 


O: Sla 
X: S2 


ISM 






070411 


O: — 
X: S2 


-0.17 + 0.06 






0:Slb 
X: S2 


Wind 







Notes. Following the scheme in Sect. l3.1.2l different criteria were applied to find a consistent description of the late-tim e optical and X-ray after- 
glow data to reveal the nature of the circumburst medium. The light-curve segments were labelled following Panaitescu| d2007L for the designated 
ID see also Table[2) ((S, J, j)=(spherical expansion, jet with sideways expansion, jet without lateral spreading), (1, 2)=(v < v c , v > v c ), (a, b)=(ISM, 
wind); see also Tabled. If there was a break in the light curve, the segments were tagged with the subscript 1 and 2. The second column lists the 
light-curve segments and all closure relations that agree with the observational data within 3<x. The third and fourth column show the correspond- 
ing measured differences in the temporal and spectral slopes, if available. Their uncertainties are lcr errors. If a criterion could not be applied or a 
solution was not found, we crossed the field out. The fifth and sixth columns show the conclusion based on all criteria, i.e., the segment label and 
the type of circumburst medium. The second last column summarises the results from the literature. 

References. (1) iRvkoff et alj d2006T> : (2) ICenko et all d2006l) ; (3) iBurrows et aT] (2006); (4) [Mangano et all d2007l) 
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GRB 


Closure relation 




Q'opt 


fix - 


fiopt 


Conclusion 


Literature 

ID 


Reference 


070802 


O: Sla/b, S2 
X: S2 


0.27 


±0.18 


0.50 


±0.05 


O: Sla 
X: S2 


ISM 


O: S2 
X: S2 


5 


07081 OA 


O: — 

X: Sla, S2 


-0.01 


±0.12 






O: Sla 
X: S2 


ISM 








O: Sib 


-0.20 


±0.05 


0.48 


±0.12 


O: Sib 




O: Sib 




0803 19B 


XiiSla, S2 
X 2 : 


Xi: S2 
X 2 : — 


Wind 


Xi: S2 
X 2 : J2 


6 


0805 14B 


0:Slb 

X: Sla/b, SZ 


-0.10 


±0.14 


0.50 


±0.13 


O: Sib 

X: S2 


Wind 


O: Sib 

X: S2 


7 


080710 


O: Sla,j2b 
X: Sla,j2a/b 
O,: — 


-0.01 


±0.09 


0.00 


±0.01 


O: Sla 
X: Sla 
Oi: Sla 


ISM 


O: Sla 
X: Sla 

1) Oi: Sla; 2) O: Sla 


8 


080721 


2 : — 
A: S2 


0.04 


±0.09 






2 : S2 
X: S2 


ISM 


1) 2 : S2 

1) A: — ; 1) A: S ia 


1)9; 2) 10 


0809 16C 


O: 
X: 

X,: S2 


-0.11 


±0.13 


0.00 


±0.34 


O: Sib 
X: Sib 

Xi: S2 


Wind 


O: Sla/b 
X: Sla/b 


11 


no i inQ a 


L): — 

X 2 : J2,j2a/b 


0.21 


±0.06 






O: jla 

X 2 : j2a 


ISM 






090102 


U: 

X: 


0.07 


±0.04 


0.02 


±0.22 


O: Sla 
X: Sla 


ISM 


(J: Sla 
X: Sla 


12 


090323 


O: Sib, jla/b 

X: Sla/b, S2,jl/2a,j2b 


-0.32 


±0.10 


0.30 


±0.18 


O: jib 
X:j2b 


Wind 


1) O: Jl, jla/b; 2) O: Sib 
1)X:— ; 2) X:Slb 


1) 13; 2) 14 


090328 


O: Sla/b, S2, jla, j2a/b 
X: Sla/b, S2,jla,j2a/b 


-0.10 


±0.09 


-0.07 


±0.21 


O: j2a 
X: j2a 


ISM 


O: Jl 
X: Jl 


13 


090726 


O: — 
X: S2 


0.37 


±0.10 






O: Sla 
X: S2 


ISM 


O: Sla 
X: S2 


15 


090902B 


O: Sla, S2 
X: Sla, S2 


0.36 


±0.04 


0.29 


±0.17 


O: Sla 
X: S2 


ISM 


O: Sla 
X: S2 


13, 14 


090926A 


O: Sla/b 
X: Sla 


-0.21 


±0.07 


0.00 


±0.08 


O: Sla 
X: Sla 


ISM 


1)0: S2; 2) O: Sla 
1)X: S2; 2) X: Sla 


1) 14: 2) 16 



Refe rences. (51 iKriihler et al.l d2008h: (61 iRacusin et ail d2008h: (71 iRossi et ail ( 120091): (81 iKriihler et al d2009h: (91 | Kann et alj 1 20101: 
(10) IStarling et alj d2009l); (111 iGreiner et alj d2009H) ; (12) lGendre et alj i20l^ K (131 iMcBreen et alj d20ld) : (14) ICenko et alj d2010l) : (15) 
ISimon et al.l OOlOl) : (16) feau et al.ld2010h 



TablelB.3l — continued 
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Appendix C: Figures of the light curve fits 



Schulze et al.: The circumburst density profile around GRB progenitors: a statistical study 




GRB 050603 

p= 1.92 ±0.24 
/3 X = 0.96 ±0.12 
Popt — 0-6 (assumed) 
-4* = 0.0 (assumed) 




Wt = 39.21 ks 
^stop — 196.5(1 ks 



4.0 
3.0 
2.0 
+? 1.0 

0.0 

" -1.0 
-2.0 
-3.0 

4.0 

"ho 3.0 
c 

2.5 
CT 0.4 

0.2 

Q. 

e o.o 

-0.2 
Jt -0.4 
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Fig. C.l. Optical and X-ray afterglow light curves of the 27 bursts that entered our sample. Upper panel: The optical data in the R c 
band are shown as dots and the X-ray data at 1.73 keV as bigger dots with an error bar in time. The light-curve fits are over-plotted. 
Upper limits are shown as downwards-pointing triangles. The grey box is the overlapping time interval of the late-time evolution. 
Vertical dotted and dashed lines indicate breaks in the optical and X-ray band. Information on the SEDs are shown in the bottom left 
(see also Table IB. lb . The given extinction, A|j°f, is the observed host-extinction in the R c band based on the deduced host extinction 



in the V-band, Ay . Additionally, we ded uced the e lectron index, p, from /3 X . The electron index is either p = 2/3 if v c < v x or 
p — 2/3 + 1 if v c > v x (e.g.. lZhang & Meszarosll2004l) . Its error was computed by propagating the uncertainty in /3 X . Middle panel: 
The flux density ratio between the optical and X-ray afterglow is shown as a solid line and its error as a dashed line for the shared 
time interval of the late-time evolution. The grey box represents the allowed parameter space of the flux density ratio (TableQ]). The 
upper boundary is the expected flux density ratio for v c < v opt , while the lower one shows the expected ratio for v c > v x . If the 
cooling break is in between the optical and the X-ray bands, the expected flux-density ratio lies be in between these boundaries. 
The expected flux density ratio depends on the electron index. Not all bursts could be corrected for host extinction. The error on 
the electron index was neither propagated into the error of the expected nor of the observed flux-density ratio. Lower panel: The 
first logarithmic derivative of the flux-density ratio, (a x - a p t ) (f), is shown as a solid curve and its error is plotted as a dashed 
line. For f/fbreak ^ 1, the first logarithmic derivative is identical to the difference in the decay slopes obtained from the light-curve 
fit (asymptotic values). Usually breaks in the light curves tend to be smooth instead of sharp. Because of this, the first logarithmic 
derivative deviates from the asymptotic value close to a break depending on the smoothness of the break. Two solid lines are plotted 
(<4 highlight the time interval when the asymptotic decay slopes were reached within lcr. The precise values are shown on the left 
and in Table[3] Within 3cr, the asymptotic difference in the decay slopes agrees either with +1/4, 0, -1/4 depending on the spectral 
and dynamical regime and the circumburst density profile. Furthermore, an envelope is drawn around expected values, +1/4, 0, 
-1/4, with a width of 0.1 to guide the eye. 
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